Novel Structural Motifs in Clusters of Dipolar Spheres: Knots, Links and Coils 
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We present the structures of putative global potential energy minima for clusters bound by the 
Stockmayer (Lennard- Jones plus point dipole) potential. A rich variety of structures is revealed as 
the cluster size and dipole strength are varied. Most remarkable are groups of closed-loop structures 
with the topology of knots and links. Despite the large number of possibilities, energetically optimal 
structures exhibit only a few such topologies. 
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Isotropic van der Waals or depletion forces favor com- 
pact packing of spherical particles, since such arrange- 
ments lead to a large number of nearest neighbor in- 
teractions. For small clusters, icosahedral motifs often 
prevailfi giving way to close packing for sufficiently large 
or bulk systems. In contrast, the energetically optimal 
arrangement of two purely dipolar particles is with the 
dipoles aligned head-to-tail, leading to chain formation. 
In systems with both isotropic and dipolar interactions, 
frustration arises from the competition between these two 
effects. 

The study of dipolar sphere systems has a long and 
controversial history, concerned mostly with the exis- 
tence and nature of their phase transitions.^ Fresh im- 
petus comes from recent experiments, in which nanopar- 
ticles of low polydispersity with a magnetic dipole have 
been synthesized and studied by electron microscopy/^ 
Dipolar particles have also been harvested from magne- 
totactic bacterial^ Colloidal systems like these provide 
great scope for exploring the influence of particle inter- 
actions on structure and dynamics both because highly 
detailed information is available from electron or confo- 
cal microscopy, and because the form of the interactions 
can be finely tuned by the experimental conditions 

Dipolar particles start to aggregate at low volume 
fractions,'^ producing chains and clusters. In this Letter, 
we identify the energetically optimal structures of clus- 
ters composed of spherical particles that interact through 
a permanent (electric or magnetic) dipole moment in ad- 
dition to an isotropic soft core and attractive tail. We 
ask how the morphology of the clusters is affected by the 
number of particles and the strength of the dipolar inter- 
action relative to the isotropic attraction due to van der 
Waals or depletion forces. 

We start with the Stockmayer potential. 
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for the total pairwise interaction energy of N dipolar par- 
ticles, where Tij is the position of particle j relative to 



particle i and p,^ is a unit vector along the dipole mo- 
ment of particle i. The units of energy and length are 
set by the Lennard- Jones (LJ) e and a parameters, while 
the dimensionless parameter /i determines the strength 
of the dipolar contribution relative to the LJ part. 

Considerable progress has been made in the field of 
global optimization in the decade since Clarke and Patey 
first searched for the lowest energy structures of small 
dipolar structures using simulated annealingii& These au- 
thors concentrated particularly on the TV = 13 cluster, 
denoted here St 13, finding a sequence of structures with 
increasing that has been verified recently.- For larger 
sizes, "intestinal" structures were reported, and the anal- 
ysis of these intricate assemblies is one purpose of the 
present contribution. We employ the well documented 
basin-hopping algorithm^ in which a Monte Carlo simu- 
lation is run on a transformed potential energy surface by 
performing a local geometry optimization at each step. 
We found that geometry optimizations converged more 
efficiently if each /ij is represented by a pair of spherical 
polar angles rather than by a Cartesian vector, since the 
latter contains a redundant third degree of freedom. 

Global optimization runs were performed for 3 < TV < 
55 particles and a grid of dipole strengths covering the 
range < /i < 6, each run initiated from a random 
configuration. Starting from the optimal structure of 
a neighboring size is often detrimental, since the global 
minimum can change abruptly with N or fi, potentially 
trapping the search in the wrong morphology. To prevent 
particles from occasionally drifting away, the cluster was 
placed inside a hard spherical container. The number and 
length of runs varied with cluster size and dipole strength 
as required to obtain reliably reproducible global min- 
ima. Since the dipolar contribution to the energy changes 
with the square of /x, the relevant energy scale increases 
rapidly as one moves away from the /i = LJ limit. 
Some care must be taken in choosing a suitable temper- 
ature kT for the basin- hopping searches, to ensure that 
the transformed surface is properly explored, while the 
global minimum structure retains significant statistical 
weight. We found that kT/e = 1 + was a useful 
starting point, though other schemes were tried in cases 
of doubt. 

Figurensummarizes the occurrence of the most promi- 
nent morphologies resulting from global optimization. At 
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FIG. 1: Structural map for global potential energy minima 
of Stockmayer clusters. Symbols denote morphologies: re- 
laxed Lennard- Jones (black filled circle), decahedral (cyan 
filled square), linear (blue cross), ring (red open circle), two 
stacked rings (purple triangle), coil (brown diamond), link 
(orange open square), knot (green plus). Structures that do 
not fall into these categories have been omitted for clarity. 



/i = 0, the well known LJ structures are reproducedji 
which, with the notable exception of the face-centred cu- 
bic 38-particle cluster, are constructed from icosahedral 
motifs. A small non-zero typically leads to a slight re- 
laxation of the LJ structure. The preference of dipoles for 
head-to-tail alignment gives rise to frustrated circuits of 
dipoles within the LJ structure, lowering the symmetry!^ 
For example, the perfect icosahedral symmetry of the 13- 
atom LJ cluster is broken by loops of dipoles that point 
around one of the three-fold axes. Even without decorat- 
ing the particles with the dipole vectors, the point group 
defined by the particle positions changes from 1^ to D^. 

At some point, the dipole contribution to the energy 
becomes large enough to favor a different structure. The 
threshold at which the change happens varies substan- 
tially with N, depending on the nature of the compet- 
ing structures. In some cases, such as 51 < iV < 54, a 
vacancy in an icosahedral shell simply finds a more fa- 
vorable site. In other cases, an entirely different class of 
compact structures takes over, such as for certain sizes 
in the range 33 < iV < 41, where there is a switch to 
a decahedral motif with at least one incomplete shell. 
Surprisingly, the dipoles circulate roughly about an axis 
perpendicular, rather than parallel, to the pseudo five- 
fold axis. The emergence of decahedra here is partly 
explained by the fact that the number of nearest neigh- 
bor pairs in decahedral structures is high, relative to the 
underlying average, for odd N close to 38,^^ making dec- 
ahedra competitive with icosahedra for alternate N close 
to this size. The exception is = 39, where the icosa- 



FIG. 2: Some multi-component global minima, (a) Stao at 
fj, — 3.6: the trivial unlink, (b) St48 at p = 3.4: two linked 
coils, 4i, and (c) Stas at fi — 3.2: the link 63. In each case, 
spherical particles are shown in the upper panel, and the un- 
derlying chain of dipoles in the lower. Components within 
each link are distinguished by color. 



hcdral structure has a stable partial second shell. This 
subtle interplay of forces is one of several examples that 
give rise to the diverse results mapped out in Fig. ^ 

In the opposite limit of high fi, the global minimum is a 
ring with the dipoles tangential to the edge for all A^ > 5. 
Compared with a linear structure, the energy required to 
bend the chain of dipoles into a ring is more than compen- 
sated by the additional contact made between the ends. 
When /i is lowered sufBciently for A^ = 10 and A^ > 12, 
the LJ contribution to the energy makes it more favorable 
to form two rings of half the size, which then stack on 
top of each other, creating one new nearest neighbor per 
particle. The threshold of /i at which this change occurs 
shows an overall increase with N, since the difference in 
bending energy per particle between a ring of A^ parti- 
cles and two rings of A^/2 particles decreases with A^. For 
even A^ the two rings are planar, forming an antiprism, 
as illustrated in Fig.|2Ia). For odd A^, the two rings differ 
in size by one particle, and an out-of-plane distortion is 
necessary to accommodate the mismatch. Since some of 
the dipole-dipole interactions are frustrated in this ar- 
rangement, the single ring becomes more favorable at a 
lower fj, for odd A^ than for neighboring even A^. The 
odd-even alternation in the boundary between one and 
two rings is clearly visible in Fig. ^ 

The antiprismatic structure in Fig. |2Ia) shows two 
clearly defined closed circuits of dipoles. While neigh- 
bors in the same chain approach somewhat more closely 
than members of different rings, making it possible to 
define chains using a distance criterion for neighbors, it 
is much less ambiguous to use an energetic definition. 
The anisotropy of the dipolar potential means that the 
propensity of a particle for "bonding" is effectively satu- 
rated by two particles with roughly aligned dipoles posi- 
tioned at its head and tail. Starting from a given particle 
i, we define the next member of the chain as the parti- 
cle J that has the lowest dipole-dipole interaction energy 
with i and is located in the half-space to which the dipole 
of i points, i.e. /tj ■ Vij > 0. The previous particle in the 



chain is likewise the particle with lowest interaction en- 
ergy that also satisfies fi^ ■ < 0. This parameter-free 
definition reliably and intuitively decomposes a structure 
into its constituent chains. In compact clusters, where 
the dipolc-dipolc interactions are not dominant enough 
for the structure to be naturally decomposed into chains, 
an attempt to define chains in this way will result in a 
collection of meaningless fragments and can be discarded 
as irrelevant. 

Analysis of structure in terms of chains allows us to 
make sense of the intermediate /i regime, where neither 
the isotropic LJ energy nor the dipolar contribution dom- 
inates. One unexpected solution to the frustration be- 
tween ring formation and condensation is illustrated in 
Fig- 13a)- In this St^s cluster, all the particles belong to 
a single closed-loop chain, but the chain has the topology 
of a non-trivial knot; it cannot be unravelled into a simple 
ring without breaking the chain. At every point in the 
knot, three portions of the chain are in close proximity, 
giving a large number of nearest neighbor interactions 
(114 pairs, based on a distance threshold of 1.35cr). At 
the same time, the structure remains open, limiting the 
bending energy, and keeping unfavorable combinations of 
dipole orientations apart. The St^s example has D2 point 
group symmetry, i.e. three distinct C2 rotation axes but 
no mirror planes. Although the projection may be unfa- 
miliar, the structure is a trefoil, the simplest non-trivial 
knot, and is written 3i in Rolfsen's notation^ meaning 
the first (in fact, only) knot whose reduced projection 
has three crossings. 

The topology of a closed-loop chain is invariant to any 
deformation (ambient isotopy) that does not break the 
chain. In contrast, the occurrence of the trefoil as a global 
potential energy minimum clearly relies on the detailed 
arrangement of the particles. Nevertheless, it is natural 
to ask how prevalent the trefoil is, and whether knots of 
different topology also arise. A reliable way to distinguish 
prime knots is through their Jones polynomials,-'^^ which 
are invariant to ambient isotopies and unique for knots of 
up to nine crossings, as well as most of ten. We find the 
Jones polynomial of a knot cluster starting from an arbi- 
trary projection of the underlying chain onto two dimen- 
sions. The so-called bracket polynomial is first derived 
by considering all combinations of splits at the crossings 
in the projection, and is then combined with a factor de- 
pending on the writhe to give the Jones polynomial. -'^'^^ 
The corresponding Rolfsen notation^" can then be looked 
up in tables 

The trefoil knot, of which Fig.|3fa) is an example, first 
appears at = 21, though it is necessarily more compact 
at these smaller sizes. This topology dominates the first 
band of knots visible in Fig.nin the range 21 < A < 38. 
However, towards the larger end of this range, a second 
topology, 5i, of greater complexity appears. Figure Ol^b) 
shows that this class of global minimum is more com- 
pact than the larger trefoils. It has three neatly stacked 
turns, with a thread through the central axis. The sig- 
nificant bending of the dipole chain is compensated by 
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FIG. 3: (a) Stag at ^ = 3.6: the trefoil knot 3i, (b) Stas at 
fi = 2.8: knot 5i, (c) St^g at fj. = 1.6: knot 8ig, (d) Stss at 
fi — 3.2: knot IO124, (e) St45 at fj. — 2.6: knot IO139, (f) St54 
at /i = 4.6: a coil with the topology of the trivial knot. In 
each case, spherical particles are shown in the upper panel, 
and the underlying chain of dipoles in the lower. The color 
changes smoothly along the chain. 



the compactness of the overall structure. 

Figure ^ shows a second band of knots, starting at 
lower fj, and rising with A^ in the range 38 < A^ < 55. The 
trefoil does not occur here, but the 5i knot is a common 
feature. We also observe two substantially more intricate 
knots whose reduced projections contain ten crossings. 
Fortunately, despite having more than nine crossings, 
neither has an ambiguous Jones polynomial,^"' allowing 
both to be identified as described above. As shown in 
Fig.EIe), one of them, IO139, resembles the 5i knot in its 
packing. The slightly wider turns now accommodate two 
threads, creating two stacks each of three turns. In the 
illustrated case of 8145, the stacks are equivalent, being 
interchanged by a C2 operation. The other 10-crossing 
knot, IO124, shown in Fig. I^Jd), more closely resembles 
the twisted wreath of the trefoil in Fig.|3Ia), but with a 
denser bundle of four turns. 

A knot with eight crossings in the reduced projection, 
819, shown in Fig. |3|[c), makes an appearance at A^ = 38 
for values of /i below those that produce the trefoil at 
the same A^. The smaller dipole moment means that 
the isotropic LJ contribution is more influential, and this 
cluster is indeed compact, having 156 nearest neighbor 
pairs, compared with 114 for the trefoil. These contacts 
come at the expense of some sharp bends in the chain of 
dipoles. 

We note that it is possible for a closed chain of dipoles 
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to possess multiple turns without passing through its own 
loops to produce a knot. Such coiled structures are also 
observed as global minima, as illustrated by the St54 clus- 
ter in Fig. |2If). Furthermore, there is the possibility of 
forming more than one closed chain in the same struc- 
ture. Such combinations are known as links, and we have 
already seen a topologically trivial example in Fig. IJl^a). 
However, non-trivial links are also encountered in the 
structural map of Fig. ^ The smallest and simplest ex- 
ample is the Hopf link, 2f , of two interlocked rings, which 
first occurs for Sti2 in the range 1.6 < fj. < 2.3 and con- 
sists of two interlocked hexagons. A less clear example is 
illustrated in Fig.[5fb) for St48, which is composed of two 
interlocking coils, each of two turns, and has overall C2 
symmetry. One three-component link topology, 63, con- 
sisting of three mutually interlocked rings, has also been 
observed, and is illustrated in Fig. |5fc). The interplay 
of factors that determine the number of components in 
a link is expected to be rather delicate, and indeed Stsa 
is a 6| link for values of fj, where St32 or St34 is a knot. 
In this case the balance may be tipped by the fact that 
the 33-particle cluster can be divided into three rings of 
equal size. 

The occurrence of such topologically exotic structures 
as global minima of the simple model Stockmayer po- 
tential was unexpected. On the other hand, given that 
knots turned out to be optimal in some cases, it is now 
remarkable that only a few topologies are observed for 
N < 55 out of the 249 possibilities with up to ten cross- 
ings. For larger clusters, where a single closed-loop chain 
would be long enough to accommodate even more cross- 
ings, the possibility of yet more complex topologies arises. 



However, the present identification by Jones polynomials 
would then be inadequate, since the polynomials are not 
unique for 10 or more crossings^ For sufficiently strong 
dipole moment, the global minimum of a large cluster 
will always be a ring to obtain the maximum number 
of head-to-tail contacts. However, as Fig. ^ shows, the 
threshold at which rings are optimal increases with size. 
In contrast, the average dipole strength at which the LJ 
structure is superseded shows no overall trend with size. 
The range of dipole strengths over which complex struc- 
tures such as knots may be found therefore widens with 
increasing cluster size. 

The fact that knots appear over a reasonable spread 
of sizes and dipole strengths suggests that there is a 
good chance of observing some of them experimentally in 
suspensions of dipolar colloids. However, while the ad- 
vantages of the Stockmayer potential include being well 
known, the LJ contribution is only one choice for the 
isotropic attraction. For some colloids, a shorter-range 
potential and stiffer repulsive core might be more appro- 
priate, and it would be interesting to see whether the 
delicate balance of nearest-neighbor and chain-like ten- 
dencies that leads to knots is preserved with respect to 
such changes. 

A number of other questions arise concerning thermo- 
dynamic stability at finite temperature, and the mech- 
anism and dynamics of self-assembly.^^ Further work is 
required to relate these properties to the overall organi- 
zation of the underlying potential energy landscaped 
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